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ABSTRACT 

With asymptotic method developed from weak turbulent theory, the 
kinetic equations for QGP are expanded in fluctuation field potential 
AJ^. Considering the second-order and third-order currents, we derive 
the nonlinear permeability tensor function from Yang-Mills field equa- 
tion, and find that the third-order current is more important in turbulent 
theory. The nonlinear permeability formulae for longitudinal color oscil- 
lations show that the non-Abelian effects are more important than the 
Abelian-like effects. To compare with other works, we give the numeri- 
cal result of the damping rate for the modes with zero wave vector. 
KEY WORDS: QGP, Kinetic theory. Weak turbulent theory Fluctu- 
ation. 
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Kinetic theory is an important method by which it is possible to investigate the non- 
Unear, non-Abehan and non-equihbrium effects in quark-gluon plasma (QGP). Some progress 
has been made in the study of the collective behaviors in QGP by the applications of classi- 
cal or semiclassical kinetic equations ||l], |2|. The medium response to external fields has been 
preliminarily discussed. However, the non-Abelian effects has not been demonstrated well 
in the works done before, only linear effects are involved in most of those works. Therefore, 
the results are similar to the Abelian case. Some authors]^ ^ have computed the non- linear 
eigenfrequency shift and non- linear Landaue damping rate with multiple time-scale method, 
but the nonlinear dielectric function, from which the properties of QGP in non-equilibrium 
state can be studied, are not discussed. In this letter, we aim at obtaining the non-linear 
permeability and the collisionless damping rate by regarding the non-linear collective effects 
as the results arisen from the turbulent fluctuations in QGP. In general, the key to calcula- 
tions beyond linear approximation is that one need to find a meaningful iterating method 
with which nonlinear equation can be effectively solved. Fortunately, a weak turbulent the- 
ory can just help us to deal with the the problem since the theory gives spontaneously the 
small quantity, the field excited by fluctuations, which is used to asymptotically expand the 
kinetic equations and Yang-Mills equation. 

In weak turbulent theory, the quantities with characteristics of collective motion can be 
decomposed into regular parts and stochastic fluctuation terms |p: 

/ = + f, f = + F, G = + G^, 

where /, /, G denote the distribution functions of quarks, antiquarks and gluons in QGP, 
respectively, is 4- vector potential, R the index of the quantities averaged over the random 
phases, namely, =< X >, and T the index of the fluctuation terms. To simplify, as is 
done in electromagnetic plasma, we let = 0, that is, the fields excited by the fluctuations 
of particle density in QGP are considered, while the background fields are ignored. 
We now outline the calculations. We start from the following equations|^, |^, 



p^D^f + ^p^d;{F^,J} = 0, (1) 

P^Dj-^p^d;{F,,,f} = 0, (2) 

p''5^G + |p^a;{F^,,G} = 0, (3) 

D^F^-=r, (4) 

f = -\j T^/^^f^f - + '^'^-fabcG,,)^ (5) 



where F^y,F^y represent the field tensors in fundamental and adjoint representations, i.e., 
F^y = F^jjTa, F^i, = F^j^Ta with the corresponding generators Tq and Tq, fahc is the structure 
constant of SU(3), Nf the number of the flavors of quarks, the 4-current density. In weak 
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turbulent approximation, is small and can be used to expand the quantities such as 
/^,Fand G^. 

At first, we deal with the case in linear approximation. By averaging Eqs(l),(2)and(3)over 
stochastic phases, we obtain the expressions for and G^, 

I^d.f^ = -igv^{[Al, f]) - |;.^{Fj„ d^f}), (6) 
P^dJ^ = -igp^{[Al, f]) + d;r}), (7) 

p^3,G^ = -zgp>^{[Al, G^]) - |p^{fJ,, ^^G^}). (8) 

Subtracting (6), (7), (8) from (1),(2),(3), respectively, we arrive at the equations for the 
fluctuation terms of the distribution functions 



-igp>^i[Aj:,r]- < [Aif^] >) - < {F^.,d;r} >), (9) 



-igpn[Aj:,F]- < [Aj:,n >) + ip>^d;{{F^,,d;n- < {F^,,d;F} >), (lo) 



p^d.G^ + igp^Aj:, G«] + Ip'^d^iP;^,, G«} 
= -igp>^i[Aj:,G^]- < [AJ^,G^] >) - lp>^d;{{F^,,d;G^}- < {F^,,d;G^} ». (11) 

In linear approximation, all fluctuation quantities are kept to the first-order in Aj^. By 
using Fourier transformation, Eqs(9),(10)and(ll) enable us to write the fluctuation terms of 
the distribution functions in terms of the regular parts of the functions and Aj^. Substituting 
the fluctuation terms into Eq(5), one gives immediately the first-order current 



s n"{k)Al. (12) 

Obviously, the covariant polarization tensor 11'^'^ should be determined in terms of 
f^:F G^. Furthermore, the 3-dimension dielectric tensor can be written as[^] 
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and eij is decomposed into longitudinal and transverse components ||8|, ^] 

e^,=eL{k)^+eT{m,-^). (14) 
Let us consider the ultra-relativistic case where po = |p|. We obtain 

..(u.,k) = 1 + ^[1 + ^1-1^1 + - (15) 

= ' + ;jw'-2k + ill - i^HH^i + - Fied"! d'') 

where = ^f^ J dpoPo^{Nf{f^ + /^) + 2NG^) is the effective mass, and G(j;) a step 
function. If the fluctuations near equilibrium are considered, Eqs(15) and (16) coincide 
with RefQ completely, and f^,f^ are determined by Feimi-Dirac distribution np, by 

Bose-Einstein distribution n^, and m = rriD = gT\J iXj: ^//^ which is the Debye screening 
mass, as determined from a 1-loop approximation in thermal QCD[|l^]. If we discuss the 
non-equilibrium QGP, the mass m will deviate from mo because of the difference of the 
distribution functions from those in equilibrium. 

In linear approximation, a possible damping mechanism for those collective modes arises 
from collisions of the plasma wave with individual particles. But one can check that the 
Landau damping has never happened as is done in Ref |]^, ^, |^] . 

We now turn to study the dielectric properties beyond the linear approximation. In this 
case, the high-order kinetic equations should be involved. Here only the second-order and 
the third-order equations are considered. With recurrent iterating method in solving the 
kinetic equations, the second- order and third-order fluctuation terms are given in terms of 
the square and cube of vl^, and G^. Furthermore, the second-order and third-order 

currents can also be determined. According to the field equation (4), one finds easily in the 
case of longitudinal oscillations that 

- \e\\^)<dg> 

f (kki)(kk2) 

= 9 J dkidk26{k - ki - ^2) n^||]^^[]^^| ^'^i^dra, u]) < abg > 

f (kik2)(kk3) 

-9^ dkidk2dkz5{k - ki - k2 - ks) , tr{Td[Ta, [n, Tc]]) < abcg > 

J |k||ki|k2||k3| 

+ / dkidk25{k - ki - k2)tr{Td[Ta,n]) < abg > LOiWk^kiM 



dkidk2dkz5{k - ki - k2 - k^)tr{Td[Ta, [n,Tc]]) 
> - < ag X be >)uJiuJ3^kMM,ka 
dkidk2dks6{k - ki - k2 - k'i)tr{Td{Ta, {n^Tc}}) 
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X (< abcg > - < ag >< be >)uJiuj3'^kMM,ks 
g'^m''^ 
16 
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X (< ahcg > - < ag X be >)uJiUJ2i^3'^kMM,k3 
4 

X {< abcg > - < ag X bc>)LJiLJ2Ld3^k,kiM,k3^ (l'^) 



dkidk2dks6{k - ki - k2 - k3){6bcSad + ^ab^cd + 5achd) 



with 



k,ki,k2,k3 



d^l vk vki vk2 1 1 
4tt |k| |ki| |k2| vk + iO+ vk2 + iO+ ' 
dO, vk vki vk2 vks 1 1 



Air |k| |ki| |k2| jksl ufc + iO+ i;(A;2 + k^) + iO+ vk3 + i0+' 

where e'(fc) denotes the hnear permeabihty for longitudinal modes, i.e., eiik) mentioned 
above, v = in u = (1, v) is the velocity of the particles in QGP, /dfi runs over all the 
directions of the unit vector v. For the sake of simplicity, we have already used the following 
marks in Eq(17) 

< ab >=< A^{ki)Al{k2) >, 

< abc >=< Al{ki)Al{k2)A^{k3) >, 

< abed >=< Al{ki)Al{k2)A^{k3)A^{k4) >, 

where <> represents averaging over stochastic phases. The average values above are called 
two-point, three-point and four-point correlations, respectively. We can approximatively 
split the four-point correlations into three pairwise products of two-point correlations in 
weak turbulent theory, whereas this procedure gives zero for the three-point correlations 
since we know < a >=< A'^{k) >= 0. These can be expressed as 

< abed >ss< ab >< ed > + < ae >< bd> + < ad >< be > 

and 

< abe >ss< a >< b >< e > + < ab X e > + < be >< a > + < ae X b >= 0. 

Here we have taken the two-point correlations as < ab >=< (ki) Aj [ki) > 6{ki + ^2) for 
a stationary uniform turbulent fluctuation |5|. Marking < A'^{k)Aj (k) > = < ab >, Eq(17) 
then can be written as 

11 f 

+ i^^'^^^^df + ^g'^fbfefade / dki < a6 > [1 + gm^ {ioT k,ki-kuk + i^l^k,kuk-ki)] 

If 11 

- -g^ / dki < ab> [m'^{—6adSbf + ^dbjedade) + ^rn'^i^adhf + ^ab^dj + ^afhd)] 

X wt^i(rfc,fci + Tfc^fci.fc < fg >= 0, (18) 
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where dabc is the symmetric constant in the generators algebra of SU(3). Since < fg >^ 0, 
its coefficient in Eq(18) should vanish. It is the dispersion equation which contains the 
contributions of non-linear currents. The total permeability is the sum of two terms 

here is the nonlinear permeability as defined by 

^df = —j dki < ab > 
1 2 2 

X [fabeffde + fnjjfbfefade{^^k,ki-ki,k+^l^k,ki,k~ki) 

— -^9^'nn''^{j^6adhf + -^dbfedade)^'^''^ k,ki~ki,k + ^k,ki,k-ki) 

- ^g'^m^gi^bf^ad + ^ab^df + <Ja/'^M)wW?(rA:,fci -fci.fc + ^kM,k-ki)]- (19) 

The above formula shows that the non-linear permeability in QGP is an 8 x 8 matrix 
in color space and contains the structure constant of SU(3), which implies the that the 
non-Abelian features have been reflected in dispersion relations. One sees that it is very 
difficult to analytically integrate (19) because of its complicated structure, the numerical 
technique is required. As an example, we only compute the eigenmodes with the same color 
for the equilibrium QGP with random thermal fluctuations, 

< ab >=< {A^{k)f > 6ab, 
< {A^{k)f >= ^[6itu - Lo{k)) + 6{lo - cu(k))]/(k), 

where /(k) characterizes the total intensity of the fluctuation oscillations with frequencies 
a;(k) and — a;(k). Therefore the nonlinear permeability matrix becomes a scalar coefficient 



-,2^2 



^ - ^ dkiI{'ki)—5{uJi-uJi{'k-i_)){u}Tk,ki-ki,k+^^i'rk,kuk-ki) 



H — xeiYTT^Lj^ ^'^ J dkili}ii)^6{uji - ui(ki))ujujj{rk,ku-kuk + ^kM,k~ki) 

5jy2 3 5 /" vr 

+ Q 2 2 ff^(~ + o) / dkiI{'ki)—6{uJi-uJi{'ki))uJUjl{Tk,ki-ki,k + Tk,kuk-ki)- 
ovr uj g 1 J uj-^ 



(20) 

The integrals in Eq(20) contain three complex factors denoted by ^j^^^q^ , 
t,(A;-fci)+iO+ • According to the Cauchy principal value formula, three delta functions 6{vk), d{vki ) 
and 6{v{k + ki)) should emerge in the imaginary parts of the integrals. However, one sees 
that the classical eigenmodes in QGP are timelike from the linear- response properties, 
so that the Chevenkov's condition [vk = 0) can not be satisfied at all. Therefore, only the 
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last factor arising from nonlinear processes makes contribution to the imaginary parts of 
the integrals. Since one can check that the condition for the nonvanishing of the imaginary 
parts mentioned above is the existence of the delta function 6{u; — uji — v(k — ki)), the 
upper-bound of the integrals over ki in Eq(20) should equal to gT approximately in the 
long-wavelength region. In this region, the dispersion relation in linear approximation is 
reduced to the form of|lC, 11 1 



and one can obtain from random thermal fluctuations [Q] 

I(k) = AttT, 

where T is the temperature of equilibrium QGP. After these considerations, one can estimate 
that the first term in Eq(20) representing pure non-Abelian effects is of order g, while the 
last two terms (Note: there are similar terms in a electric-magnetic plasma) are of orders g^ 
and g^, respectively. This implies that the pure non-Abelian effect is more important in the 



nonlinear response theory. To compare with other works [12, |T^, as an example, we prepare 
to work out the result for k = modes. Calculating of in Eq(20) to the order in g, we 
have 



a;,k) = ^^fpl ^dco,d\k^\I{k^)^6{L0,-u;^{k^)) 
lovr^LiJ'^ J ujf 



i I -^1 |ki 

- TT^ea-^)]. (21) 

ojwjlkil |ki| 

We immediately obtain the numerical result of the imaginary part of 

Ime^ -0.72g. (22) 
Then the damping rate 

7 = 



del 

^2r 



= 0.25g^T. (23) 
For pure gluon gas, we have 

7^ = 0.20g'^T. (24) 
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Eq(24) coincides with the results based on the resummation of the hard thermal loops [p!2| 
and the asymptotic calculations with the multiple time-scale method fl^ . In fact, more sit- 
uations can be discussed from our method as long as one knows the knowledge of the mean 
distribution functions denoted by index i? in a nonequilibrium plasma. 

Finally, we summarize our paper. 

(1) It is an effective method to apply the weak turbulent theory to the study of non- 
linear problems in QGP. In principle, we can derive any permeability and damping rate 
(collisionless) beyond linear approximation. In this paper. We compute only the effects 
coming from the second-order and third-order currents in the processes of three or four 
waves. The formulae of the nonlinear permeability and damping rate are obtained. For 
k = modes, the numerical result of the damping rate coincides with the results obtained 
from other methods. Our results have already show that the contributions of pure non- 
Abelian terms to nonlinear effects are more important than those of Abelian-like terms. 
Although we know that the resonant condition of eigenwave with individual particle, i.e., 
vk = 0, is not satisfied, the resonance of the particle with the secondary wave coming from 
coupling of eigenwaves in nonlinear processes is possible, i. e. the condition, vlS.k = 0, can 
be satisfied. 

(2) Three wave correlations make no contribution to the nonlinear effects in our results 
because of the special treatment in this paper. In fact three wave processes must be more 
carefully expanded with four wave process This problem will be discussed in detail in 
another paper. 

(3) It is emphasized that weak turbulent theory not only provides a method to study 
the nonlinear problems in QGP, but also can advance such studies into an important and 
specific physics field. 

This work is supported by the National Nature Science Fund of China. 
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